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Recent theoretical and experimental developments in the field of electron vortex beam physics 
have raised questions on what exactly this novelty in the field of electron microscopy (and other 
fields, such as particle physics) really provides. An important part in the answer to those questions 
lies in scattering theory. The present investigation explores various aspects of inelastic quantum 
scattering theory for cylindrically symmetric beams with orbital angular momentum. The model 
system of Coulomb scattering on a hydrogen atom provides the setting to address various open 
questions: How is momentum transferred? Do vortex beams selectively excite atoms, and how can 
one employ vortex beams to detect magnetic transitions? The analytical approach presented here, 
provides answers to these questions. 0AM transfer is possible, but not through selective excitation; 
rather, by pre- and post-selection one can filter out the relevant contributions to a specific signal. 


I. INTRODUCTION 

The field of phase vortices in beams has a relatively 
short but vibrant history, beginning in optics and end¬ 
ing in all forms of wave mechanical systems. A good 
overview of theory and applications in optics can be 
found in refs. [I, 2]. Interesting aspects include the quan¬ 
tization of orbital angular momentum (0AM) [3], the 
opto-mechanical effects of these beams [3, 4], quantum 
information, computing, and communication uses [5, 6], 
and of course scattering applications [7-9] . Various ways 
of obtaining electron vortex beams have been demon¬ 
strated [10-17] since their inception on paper in 2007 [18] 
and experimentally in 2010 [12, 19]. Several investiga¬ 
tions into the scattering aspects of electron vortex beams 
have been published [20, 21], and this work aims to add 
to that list a thorough and general analysis of inelastic 
electron vortex scattering. 

The main questions we investigate in this paper all per¬ 
tain to the detection of and mechanism of 0AM transfer 
of a fast vortex electron to an atomic system. The results 
presented in this paper describe the result of a detailed 
investigation into the quantum dynamics of 0AM trans¬ 
fer from an incoming electron vortex beam to an atomic 
electron. It is an attempt to describe what the incoming 
and outgoing 0AM and transverse momentum do to the 
scattering amplitude, and how these can be exploited in 
discovering what happened during a scattering event. In 
what follows, inelastic scattering of fast electron Bessel 
beams is calculated for the model system of a hydrogen 
atom, of which several transition amplitudes are explic¬ 
itly calculated. Although we consider a spherically sym¬ 
metric atomic system, we prefer to discuss everything 
relevant to the beam axis, and we treat transitions with 
0AM transfer in that view. Special detail will be given to 
how and when selection rules are present, and if they are 
not, we describe the physical reasons behind their vanish¬ 
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ing outside of special circumstances. We put the empha¬ 
sis on what vortex beams contribute to the scattering, 
and what they do not. To this effect, two alternative ap¬ 
proaches are pursued, allowing for a sufficiently nuanced 
and accurate report on the properties of electron vortex 
beam atomic excitations. Critical are the concepts of pre- 
and post-selection on the scattered wave, and special at¬ 
tention is given to transitions in which the beam’s 0AM 
is transfered to the atomic electron. Effects related to 
relativity, spin or recoil of the nucleus are ignored. This 
treatment is also more simplistic than the more advanced 
treatments based on the density matrix [22, 23], and re¬ 
late more closely to the earlier theoretical advances in 
inelastic electron-atom scattering [24-26] . 

The next section will quickly review various inelastic 
scattering concepts, and puts special emphasis on the 
details often omitted in a textbook treatment of quantum 
scattering on atomic systems which are important here. 
The third section presents the two methods developed 
to attack this problem analytically. The fourth section 
is a discussion of the most important results. The final 
section presents our conclusions. 


Notation 

In order to remove clutter from crucial steps in the cal¬ 
culation, strict notational conventions are defined here. 
Real space coordinates are denoted r = (cc, y, z) = 
{r±,Lp,z). Basis vectors are written as where i is the 
relevant coordinate. Momentum space coordinates are 
denoted k = {kx,ky,kz) = {k±,(j),kz). It is often ad¬ 
vantageous to split any 3D vector v into a z-component 
Vz and a perpendicular component i>_l = {vx,Vy) with 
size uj_. Primed variables are denoted r' = {x',y',z') = 
(r]_, (/?', z') and similar for k' . Partial derivatives are de¬ 
noted dx = -^ and are taken to act on everything to the 
right of the symbol. Primed coordinates and quantum 
numbers refer to the final (outgoing) state. Momentum 
transfer is denoted q and is equal to k — k' . 0AM transfer 
is denoted Am and is equal to to — mb 
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II. A CYLINDRICAL VIEW ON INELASTIC 
ELECTRON PLANE WAVE SCATTERING 


A. Inelastic scattering of an electron on an atom 


In this section, various entities are defined and plane 
wave scattering theory is quickly reviewed. Certain often 
forgotten aspects which are important to the cylindrical 
wave case are brought forward. 


1. Interaction potential 


The interaction potential for an electron scattering on 
a single-electron atomic system is given by 


P = e2 




( 1 ) 


Here, r' is the lab coordinate for the scattering electron 
with respect to the nucleus, and r is the atomic electron’s 
coordinate relative to the nucleus (which is located at the 
origin). Extending this to multiple single-electron atomic 
states is formally trivial: summing over every electron’s 
coordinate suffices, so that r —>■ The relative dis¬ 

tance \r — r'\ will complicate the cylindrical scattering 
treatment in the next sections. 


2. Single electron atomic wave functions 

To make the role of the atomic system’s 0AM appar¬ 
ent, we assume that the (projected) 0AM is a good quan¬ 
tum number, and separate the atomic wave functions into 
an azimuthal part and the remaining {r±, z) part: 

Am'if 

For a hydrogen-like atom with (effective) nuclear charge 
Z and atomic radius = aome/n (p. is the reduced mass 
of the atomic system), the orthonormal wave functions 
have the following form: 




( 3 ) 

Where r] = 2Z/ (na^), and n = 1, 2, 3 . .., / = 0,... n — 1, 
and m = are the radial, angular and magnetic 

quantum numbers, respectively. The 0AM of a state, 
determined by the operator L, is given by /, and the pro¬ 
jected 0AM, given by and the one of primary interest 
here, is given by m. In what follows, we shall assume 
/i Ri rrie and so r; oq, but keep Z so that any depen¬ 
dence on atomic number is clear. 


3. Scattering amplitude 

The inelastic scattering amplitude in the first Born ap¬ 
proximation can be written as follows: 

fMk',<f] = --^{k'\{f\Vm), (4) 

AttTi 

where N takes care of the normalization of the scatter¬ 
ing electron’s states [27], \i) and [/) are the initial and 
final atomic states, |$) is the incoming state (tradition¬ 
ally taken to be a plane wave). The outgoing plane wave 
is defined as: 


{r\k') 


eifc'-r- 

(27r)3/2‘ 


( 5 ) 


The scattering amplitude determines the modulation of 
the outgoing (scattered) spherical wave so that the total 
wave function obeys the following relation: 


'i’ir) = <^{r) + ff,[k',<^]^—. (6) 

The scattering amplitude, Eq. (4), depends on the atomic 
states involved in a certain transition, and thus the prob¬ 
ability of scattering an electron in a certain direction de¬ 
pends on these states. 

Filling in Eq. (4) with Eq. (1), we obtain: 




rUee^N 

2Trh^ 


-Z(/|*)(fc'|i|$)) . 


( 7 ) 


4- The final state 

The outgoing electron’s momentum, k' is related to 
the incoming electron’s momentum k through the energy 
transferred to the atomic system, AE, as follows: 

^,2 ^ ^2 _ ( 8 ) 

The scattering angle B can be introduced by substituting: 

= 1^ E k'^ — 2kk' cos 9. (9) 

Using the above with Eq. (7) and Eq. (3), one can 
plot the angular dependence of the scattering amplitudes 
for various hydrogen transitions. These are shown in 
Fig. la. The different final states (s/p/d/...) result in 
specific angular regions in which some dominate or are 
suppressed with respect to others, which implies some 
transitions can be filtered out roughly by post-selection 
on 9 (which directly corresponds to outgoing transverse 
momentum). This can be exploited to map anisotropic 
bonding in crystals [28, 29]. 
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a) Plane wave scattering amplitudes 
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b) Is —1 2p cross sections 
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FIG. 1. (color online) a) Plane wave scattering amplitudes 
for 120 keV electrons for various transitions from the ground 
state. One can see the relative size of the amplitudes shifts 
depending on the exact scattering angle in question, b) Sep¬ 
arate cross sections for 2pz and 2p± excitation. The typical 
Lorentz profile appears [28]. The total Is —I 2p cross section 
is the sum of the other three. Note the lower 0 scale is a 
factor 100 times smaller than the one in the top graph. This 
tiny scale is due to the relatively small AE for pure hydrogen. 
More realistic systems will show the same features at a much 
larger scale. 




We will now consider a ground state excitation to a 
state with fixed (n', I'), e.g. the 2p orbitals, and see if the 
outgoing wave’s 0AM can be used to provide more details 
of the final state. The question we ask ourselves is: what 
determines scattering to a particular (projected) 0AM? 
This example final state consists of three substates: 2pz, 
2p+, and 2p_. The orientation of these final states is de¬ 
termined by the direction of the z' axis, along which their 
angular momentum is projected to give quantum num¬ 
ber m!. This is crucial to the whole analysis. Using the 
Fourier translation theorem, the plane wave scattering 



FIG. 2. For plane wave scattering, the orientation of the final 
state depends on the scattering angle 9 and the transition 
energy AE as shown in Eq. (13). The figure also shows the 
azimuthal angle, (j)' of the outgoing wave with respect to a 
fixed set of axes. 


amplitude becomes 


9rr? piq r' _ v 


( 10 ) 


This can be calculated in position space by aligning the 
z' axis along the momentum transfer vector q\ 


//^(q) = - 


2mee^ 


J d?r' 


J^qr' cos 9' _ ^ 




( 11 ) 

The final states are thus automatically quantized along 
the direction of q, when expressed in the r' coordinate 
system. [30] Furthermore, due to the integration over the 
azimuthal coordinate p' in Eq. (11), only the tilted states 
for which Am = 0 contribute. In our example of ground 
state excitation, this means that for a certain 
only the state with m' = 0 gives a non-zero contribution. 
Indeed, the transition amplitude for Is —> 2p is given by 
following expression: 


x,4>' 


{‘^Pz 


^|ls)= (- 

r \ao 


12iV2 


(—Y 

\ 2ao J 


3 • 


( 12 ) 


The subscript x expresses that the quantization axis is 
tilted over an angle of y, defined by 


q± k' sin 9 


(13) 


and rotated around the beam’s z axis by an angle (f)'. This 
scattering geometry is shown in Fig. 2. The reader can 
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verify that Eq. (12) is zero for the other tilted substates 
(in our example, these are 2p+ and 2p_). This tilted state 
(z' is rotated with respect to the beam’s z axis) can be 
projected onto untilted states which are quantized along 
the beam direction. 

As the atomic orbitals are determined by a radial func¬ 
tion multiplied with an orbital angular momentum eigen¬ 
state \l,m), their rotation properties are fully determined 
by the Wigner D-matrix: 

{l,m\R{a,j3,-t)\l' ,m') = /3, 7 ). (14) 

Rotating an angular momentum eigenstate over the an¬ 
gles a, /3, and 7 (in the z-y-z convention) turns it into a 
sum of states of equal I, but all m: 

-\-l 

R{a,P,l)\l,m') = ^ i:»®^,(a,/3,7)|Z,m). (15) 

m— — l 


For the rotation of an m' = 0 state oriented along the 
vector q as in Fig. 2, the matrix elements take on the 
following form: 


^U(<^'^X,0) = e-™'^'dL.o(x) 


= e 


{I — to )! 
{I + to )! 


^z'"(cosx), 


(16) 


where P;"* are the associated Legendre polynomials. 

As an example, consider the Is —>■ 2p transitions of 
the hydrogen atom, which can readily be calculated by 
filling in Eq. (11) with Eq. (3). With respect to the beam 
direction, the excited 2pz state is tilted by the angle % 
as shown in fig. 2. Applying Eqs. (15) and (16), for this 
example, one obtains: 


\‘2Pz)x,<l>' 


Here, 


(e’‘^'| 2 p_) -e ‘'^'| 2 p+)) -hcosx| 2 pz). 

(17a) 


sinx = 


5 

q 


cosx=—, (17b) 

q 


This expression gives the projection of the tilted state 
| 2 pz)x, 0 ' into its untilted components (quantized with re¬ 
spect to the beam’s z direction). 

The explicit form of the projection of the q-oriented 
2^2 state onto a beam axis-oriented 2p state, e.g. 
{2p±\2pz)x,<i>', in Eq- (17), allows us to write the scat¬ 
tering amplitude to e.g. the 2p± states directly using 
Eq. (12): 


fPW 

J ls—¥2p^ 


2 TOee^ 


2mee^ 



(2p±| —|ls) 


{2p±\2pz)x,4>'^ ^2 

Ti<^' -12iTOee^qL 




q^ + 


\ 2 ao J 


1 3 ■ 


I Is) 
(18) 


This is obtained by inserting a complete I = 1 basis and 
remarking that the transition element is only non-zero 
for \2pz)x,4:', as discussed above. Note that the out¬ 
going wave has acquired the opposite 0AM of the ex¬ 
cited atomic state, and thus 0AM has been transferred 
to the bound electron. The above result is not new, but 
it usually appears as a cross section. More specific forms 
are exploited in e.g. momentum resolved EELS experi¬ 
ments [31, 32], where dynamical diffraction also plays a 
large role in the final distribution of scattered electrons. 
The 2pz excitation, where no 0AM is transferred, can be 
calculated directly as well: 


fPW 
J ls—^2pz 


— 12^/2imee^ qz 





(19) 


Both Eq. (18) and (19) are shown in Fig. lb. These scat¬ 
tering amplitudes can be separately observed if one filters 
the outgoing wave on its 0AM. Techniques to achieve this 
for electromagnetic waves exist [33, 34], but in electron 
optics these are still in development [35, 36]. 


B. The vortex beam basis state: Bessel beams 

The simplest form of vortex beams is provided by the 
solution of the Schrodinger equation in cylindrical coor¬ 
dinates: 


„U(p pifciZ 

{r\kj) = il^kAf) = 

This exact solution encompasses the beam features that 
interest us: the quantized (projected) 0AM M, and the 
longitudinal and transverse momenta Hkz and hk±. One 
can define an opening angle as a = ikAkz) which 

is the angle the momentum vector makes with the z axis. 
The energy of this field-free state is independent of its 
0AM: E = A{k^ + kf)/(2,m). The Bessel beam can be 
written in terms of its momentum components, which 
shows that this state is a ring of tilted plane waves in 
momentum space: 

Mr) = {-iY I r ( 21 ) 

This representation was used to calculate the elastic 
Coulomb scattering amplitude in an earlier work [27], of 
which the results will be useful here in simplifying some 
of the equations. 

Bessel beams are basis states, much like plane waves. 
It is impossible to create true Bessel beams, even in a lab 
setting. Only approximations can be realized, which ad¬ 
mittedly show the various properties of real Bessel beams 
in a limited way [37]. Nonetheless, they provide a good 
basis to calculate scattering amplitudes because they en¬ 
code features such as cylindrical symmetry, convergence 
angle (transverse momentum), and 0AM in a natural 
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way. Note that for a convergent beam, the wave function 
of the probe in real space is nothing more than a coher¬ 
ent superposition of Bessel beam basis states, integrated 
over the aperture radius. For a non-vortex probe, this 
reads: 

POO 

4'(r) oc / dk±A{k±)k±Jo{k±r±). (22) 

Jo 

Here, A{k±) describes the aperture (for example as a step 
function). This is nothing new, and often called the Han- 
kel transform. What this means is that the scattering 
amplitude of a real convergent beam can be obtained 
from that of a Bessel beam by integrating over a k± range 
given by the aperture radius. For a vortex probe, one can 
write: 

nOO 

^^(r) oc e‘^‘^ / dk^A{kj_)k±Mk_rr±), (23) 

Jo 

which is nothing more than the ^-th order Hankel trans¬ 
form [38]. 

III. INELASTIC VORTEX BEAM SCATTERING 
A. General formulation 

We consider a perfectly centered vortex, so that the 
atom lies exactly on the beam’s 0AM axis. Beam dis¬ 
placement can be taken into account by using the Bessel 
addition formula (see Eq. (25)) for the incoming beam 
in the transverse plane. Any non-zero displacement of 
the beam (which is to be expected in any realistic situa¬ 
tion) will thus introduce a progressively larger number of 
0AM modes significantly contributing to the scattering. 
The extra modes contribute with a factor determined by 
the transverse momentum, displacement, and the 0AM 
of that contribution through Je-fj,{k±rQL), where I is the 
pure 0AM mode displaced over a distance tol, and \i is 
the 0AM of a mode introduced by the displacement. For 
roL = 0, all other contributions disappear. These addi¬ 
tional 0AM modes will coherently contribute to the final 
scattering amplitude, and thus interfere upon calculating 
the differential cross section. 

Following up on ref. [27], one can calculate inelastic 
scattering amplitudes by replacing jd)) in Eq. (7) with a 
Bessel beam |fc,.^), Eq. (20). Ignoring the (purely elastic) 
term (which was calculated in ref. [27] and repeated in 
App. A), this leads to: 


y 



FIG. 3. In-plane geometry as used in Eq. (25) for the shift to 
the atomic electron’s coordinate. 


integral as to not involve the atomic wave functions al¬ 
ready. The relative distance \r — r'\ is difficult to sub¬ 
stitute directly in this expression, due to the presence of 
the transverse coordinate. There are two ways to pro¬ 
ceed: use the Bessel addition theorem to mathematically 
displace the Bessel beam to this relative coordinate, or 
introduce the Eourier representation of the Bessel beam, 
Eq. (21). Both methods lead to interesting physical in¬ 
sights, and thus both are treated below in detail. 

Note that much like the elastic results presented in a 
previous article [27], any of these scattering amplitudes 
can be summed/integrated over a certain range of k±, k^ 
and/or 0 to more closely resemble real wave packets. This 
allows for modelling of a focused electron probe with a 
specific convergence angle (including annular apertures), 
or a certain collector angle range for comparison with the 
experiment. 


B. Displaced Bessel beam representation 

The relative coordinate problem described above can 
be worked around by applying the Bessel addition theo¬ 
rem, with which a displaced Bessel beam can be replaced 
with a series of weighted centred Bessel functions over all 
orders. In our case, the theorem states [39]: 


X f dVe-''“'-'^—e’Vj,(fc 

J \r-r'\ 

(24) 

The relative atomic coordinate of the atomic electron, 
= r — r', plays a central role in how to solve this 
problem. We would like to substitute this in the inner 


pifczz' pi^v' 

+ 00 


= e 


ikz 


Jfj.{k±r_L)e'^ 




fi——oo 








(25) 

where the relative coordinate is denoted = r — r'. 
This is illustrated in Fig. 3. 
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FIG. 4. Geometry for the outgoing plane wave’s Jacobi-Anger 
expansion in Eq. (27). Here, 7 is the angle between —k± and 

r±. 


Applying this to Eq. (24), one obtains: 

X - fl) 

+ 00 

+ 00 

= - E fel{k';k,e-^^)Mf.{k',k,^l) 

fi— — oo 

(26) 

In the above expressions, is the Coulomb scattering 
amplitude as given in App. A. This form resembles the 
plane wave result quite closely: a Coulomb scattering 
factor ( 1 / 9 ^ for plane wave incidence, assuming the form 
of /gY here) multiplied by a matrix element A4/i over the 
two bound states involved. This matrix element will take 
on a simpler form below, although the functional form of 
the “matrix element operator” is not the same as in the 
plane wave case, i.e. it is not simply e“^'’”. 

The outgoing plane wave can be written as a series 
of Bessel functions using the Jacobi-Anger theorem (see 
also Fig. 4): 

+ 00 

A= —00 

which makes the remaining hidden azimuthal angular de¬ 
pendence explicit, exposing the possibility of selection 
rules on (projected) 0AM quantum numbers. The ma¬ 
trix element then becomes: 

- 1-00 ^ 

( 28 ) 


fj,{k';k,i) = - 


mee^N 

2Trh^ 


+ 00 

E 

fl— — CC} 


Assuming the initial and final states are eigenstates of 
the projected 0AM operator {i.e. they have the form 
of Eq. (2)), the azimuthal factor can be integrated out: 


+ 00 


Mfi = E (l^\Mk±r±)J\{k[rJe^‘‘^^\a) 


X— — C 


I 


X dtp- 




2tt 


= {l3\Jti{k±rj)J^+Am{k\rj_)e^‘>^"\a) 


The change in the atomic electron’s (projected) 0AM, 
Am = m — m', is the 0AM transferred by the scattering 
electron in the collision. The reduced matrix element is 
equal to: 

MfSa = (/3|'/^(fcj_?jJ-^M+Am(fclrj_)e*'^"^|a). (30) 


This is now the only unknown. Note that due to the form 
of (see App. A) and Mfi (see above), the scattering 
amplitude’s {i.e. the outgoing wave’s) azimuthal depen¬ 
dence is exactly ^ which implies the transfer of 

the scattering electron’s 0AM to the atomic state. 

For the hydrogen wave functions (oriented with respect 
to the beam axis, and not the 0-dependent momentum 
transfer q as in Sec. IIA 4), Eq. (30) becomes a triple 
Bessel integral which can be formally solved using a col¬ 
lection of tricks. The result is unfortunately unwieldy, 
and gives no further insight into the physics of the prob¬ 
lem. 

The next section provides a alternative analytical 
treatment that can be used to calculate the scattering 
amplitude for any specific transition, both numerically, 
and analytically. This result is still useful though when 
analysing central scattering, for i.e. 6 = 0, Eq. (30) 
simplifies significantly, and due to the generality of the 
obtained expression, an 0AM reciprocity theorem can be 
deduced. These cases are discussed in Sec. IV B and IV A. 


C. Fourier representation 

The Fourier representation of a Bessel beam given by 
Eq. (21) provides a solution for the problem of the rel¬ 
ative coordinate described in Sec. Ill A. Starting from 
Eq. (7) (with N = (27r)®/^ as in ref. [27]), and using 
Eqs. (1) and (21), one arrives at the integral in question 
{£' is the outgoing beam’s 0AM, obscured by and 
the integral over the azimuthal Fourier coordinate): 

fj^{k',£';k,£) = ^ I d0e‘%Y'^(g), 

=-E(-i/ / d^—(/|e“^--Z|*). 

nh J Q 

(31) 

where the plane wave scattering amplitude, fj^{q) is 
given by Eq. (11). This can be analytically calculated 
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FIG. 5. (color online) For cylindrical waves, and by extension 
vortex waves, there is no single momentum transfer. Instead, 
the various plane wave components of the cylindrical wave 
contribute to the final outgoing wave vector component. This 
is made explicit by the integral over the various momentum 
component vectors in Eq. (31). A Bessel beam of ^ = 1 is 
shown, and the relative phase of the various contributions is 
shown by the hue of the arrows. 


FIG. 6. (color online) Top-down view of Fig. 5, showing the 
transverse plane scattering kinematics for an ^ = 1 Bessel 
beam. 


using the contour integration technique first described in 
ref. [27]. Several descriptive examples are treated explic¬ 
itly below. 

Eq. (31) has a direct physical interpretation, which is 
illustrated in Fig. 5. For a certain scattering angle 9 
(with respect to the beam’s principal direction) and en¬ 
ergy E' = /{2me), each component plane wave of 

the incoming Bessel beam will coherently interfere and 
contribute to the outgoing wave. This makes it impossi¬ 
ble to define a single momentum transfer q and thus also 
a unique final state orientation as discussed in Sec. IIA 4. 
So, in general, a (()-dependent projection is required if one 
wants to discriminate final states defined with respect to 
the beam’s axis. 

To calculate the Bessel beam scattering amplitude an¬ 
alytically from Eq. (31), one must first obtain Eq. (11), 
and then apply the following substitutions: 


2 = 

COs((/> — (j)') = 


+ k'^ + — 2Aqfc]_cos {(j) — (j)'), (32a) 

(32b) 


tan 


1 -z^ 

1 +^’ 


(32c) 


One can then extend the integration over z (over the 
real axis) to an appropriate contour (such as an infinite 
semi-circle as previously used for elastic Coulomb scat¬ 
tering [27]), and unleash the residue theorem. As such, 
the Is —>■ Is (elastic) scattering amplitude is: 


2m^e^ (R 2 - Ri\ 

= -—(-■)* 

alRlRl + 2^2 ^^2 + i?2 ^ 4l£li?ii?2) 

^ a^RlRl (33) 

+ (^J-~ ^ 1 )^ 

-^2 = + ^ 1 )^ 


d{(j) - (j)') 



2 dz 


1 -I- 2:2 ■ 


(32d) 

(32e) 


Note that the outgoing 0AM, £', is equal to the incom¬ 
ing 0AM because none was transferred to the atomic 
electron. Inelastic scattering amplitudes can also be cal- 
















culated; below is the result for the Is —>■ 2s transition: 


__y2^ /Zy , ^ 

[aoj ^ ^ {R 1 +R 2 ) 

^ 3{Rf + R*) + 6\i\RiR2{Rl + Rj) + 2(1 + 2e^)RlRj 

2a^) + + (fcj- - ^1)^ 

(34) 

For non-spherical transitions (i.e. where Am ^ 0), 
care must be taken about the final state orientation as 
discussed in Sec. IIA 4. Because there is now no unique 
momentum transfer q, and there is both incoming and 
outgoing transverse momentum, the rotation over (j)' be¬ 
comes one over (jiq as shown in Figs. 5 and 6. As an 
example, we treat the Is —^ 2p transitions, for which 
the plane wave scattering amplitude is given by Eq. (12). 
For the rotation over ^q, one can employ the complex 
representation of vector addition: 

- k'J^' 

= e*-^' . (35) 

This, together with Eq. (17b), makes clear that the in¬ 
tegrand in Eq. (31) is only dependent on (j) — (j)' and q 
(which itself is only a function of that same variable, see 
Eq. (32a)). Eor the Is —>■ 2p± transition, this gives: 


fV 
J Is— 


6imee^(-l)^ f ^ 
nh^ \ao) 


X 


J d(</. - (j)') 





(36) 

Note that the outgoing wave has lost/gained 0AM: the 
beam electron transfers 0AM to the atomic state. The 
integrated form is quite lengthy, so it is not shown here. 
The remaining 2pz state has the following scattering am¬ 
plitude: 


fu^2p^ — 


6i-\/2mee^(—1)^ 
■nh^ 



X y d((/> - </)') 


1 



(37) 

The analytical expressions for these scattering ampli¬ 
tudes can be explicitly written down, but are too complex 
to present here and don’t give much physical insight. 

As an alternative to analytical calculation, Eq. (31) 
can be calculated numerically as a function of scatter¬ 
ing angle by using Eqs. (8) and (17b), and the following 



FIG. 7. (color online) Representation of the symmetry ex¬ 
posed in Sec. IV A between incoming and outgoing vortex 
waves in the case of incoming and outgoing plane waves, re¬ 
spectively. A specific transition has the same scattering am¬ 
plitude at the scattering angle 6 corresponding to the inverted 
geometry as shown here. 


substitutions: 

q^ = k‘^ + k'^ sin^ 0 — 2k±k' sin 6 cos {4> — (j)'), (38a) 

ql = {kz — fc'cos0)^, (38b) 

and integrating over d{(j) — (j)') numerically. The func¬ 
tions in question are all in all well-behaved and standard 
numerical integration methods should have no issues with 
them. The option of simple numerical integration is very 
useful if one wants to explore a different basis set of state 
wave functions, as long as they can be quantized in 0AM 
as in Eq. (2), and their scattering amplitude can be writ¬ 
ten as a function oi (j) — 4>' explicitly. 

Figures 8, 9, and 10 show the scattering amplitudes 
for several characteristic transitions and different input 
beams. In the limit £ = 0, k± = 0, which is the red, long- 
dashed line on each of the plots, the results of Fig. lb 
are recovered. The cross sections are discussed in more 
detail in Sec. IV C 


IV. DISCUSSION 
A. OAM Reciprocity 

In the cylindrical scattering amplitude, Eq. (26), one 
can consider two specific cases: incoming and outgoing 
non-vortex states. Eor a fixed transition with a certain 
Am, an incoming plane wave (so with £ = 0, k± = 0) 
will result in a scattered vortex wave with I' = Am. Ad¬ 
ditionally, the summation in Eq. (26) can be performed 

using Jfj.{k±r±) Using the explicit form of 
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(see App. A), one obtains the relatively simple expres¬ 
sion: 




Am A^m 


Ze^ 


kl 


(39) 

Doing the same, but taking the outgoing wave to be a 
plane wave {£' = 0, fci = 0), one obtains a very similar 
expression; 


fj^{e = 0-k,i)= (40) 

X (,8|JAm(fc±n.)e''^"''|Q;). 

What these two equations say is quite simple: an inci¬ 
dent plane wave can gain/lose 0AM by exciting a tran¬ 
sition, and will then be scattered to a certain angle 
9 = tan“^(fc/_/A:'). However, a vortex wave with trans¬ 
verse fee will have the exact same probability of excit¬ 
ing that transition when scattering to 0 = 0. Note also 
that only when an incoming vortex beam’s 0AM matches 
the transition’s Am, can it be scattered to 0 = 0 (see 
Sec. IV B). This can be seen as a form of reciprocity [40- 
42], but in this case for inelastic scattering including 
0AM. 


B. Central scattering amplitudes 


Instead of considering the full 0-dependent scattering 
amplitude, one can also only consider the 0 = 0 case, for 
which the analytical expressions are much more tractable 
than for the complete calculation. Eq. (40) can be calcu¬ 
lated directly using known integrals [43] of a set of Bessel 
K and J functions. Alternatively, the analytical method 
from Sec. Ill C can be used to the same effect. 

For Is —>• Is, this gives: 


/i;^i .(0 = O;fe,£) 


2aoe'^meStfi 



/ 

1 + 
V 



(41) 


Note the extra term which comes from the first term 
in Eq. (7). If one compares this with the plane wave 
result in ref. [44], one immediately sees the symmetry 
of this expression if one replaces k± with fc]_ and realizes 
q^ = k'l + q^ in this situation. This result can also be 
obtained by setting = 0 in Eq. (33). 

For the Is —?► 2s transition, the following result is ob¬ 
tained: 


ru^,,i0 = O;k,i) 
8-\/2mee^ / Z 
V^o, 




kl 


\ 2ao J 


(42) 


The Kronecker delta ensures the outgoing wave does not 
gain or lose 0AM with respect to the incoming one, as 
only the i = 0 mode can be non-zero at 0 = 0. 

By only considering central scattering, one can se¬ 
lectively measure a specific transition by preselecting a 
proper incoming vortex state. This can be shown by con¬ 
sidering an incoming £ = ±1 vortex beam and Is —^ 2p± 
atomic transitions. The relevant central scattering am¬ 
plitude is given by: 


rV 

J ls—¥2p± 


(0 = 0;k,e) 


12imee^ 5t^±i 

kl + ql 



k± 



(43) 

For comparison, the central scattering amplitude for the 
Is —>■ 2pz transition is given below: 


fYs^2p^ — Oj k, £) 

12\l2mee‘^ bi 


kl 


qz 


kl 


Y 2ao J 


3 ■ 


(44) 

Note how the roles of k± and qz are reversed with respect 
to Eq. (43) (typical behavior for these p-character final 
states), and that the central Pz scattering cross section 
is twice as large for the same parameters. Take special 
note of the strict selection rule for 0 = 0 expressed by 
the kronecker d’s, showing that the on-axis intensity for 
these transitions will be non-zero only for the right in¬ 
coming beam. Due to the summation in the full expres¬ 
sion, Eq. (26), a mixture of outgoing vortex waves will 
generally be emitted from each scattering event regard¬ 
less its Am. The rate at which these various components 
contribute is determined by the weighting expressed by 
that equation, which is not trivial. 

Lastly, it is important to note that the three 2p states 
considered here are degenerate unless a magnetic field is 
applied, which adds a Zeemann energy, splitting the non¬ 
zero 0AM levels with the magnetic field. Larger fields 
will also induce spin-orbit coupling in an atomic system, 
further complicating the wave functions and interactions 
involved. [45, 46] 


C. Hydrogen scattering amplitudes 

Several properties of the scattering amplitudes pre¬ 
sented in IIIC, Eqs. (33)-(37), are immediately apparent 
by only looking at the equations themselves: 

• Higher energy levels {i.e. larger n) introduce higher 
order i dependence. 

• Transitions between spherically symmetric states 
(for which the scattering amplitude is independent 
of 4>) only depend on the magnitude of i, and not 




















10 


K| = 0 




d [mrad] 

FIG. 8. (color online) Is —>■ Is (elastic) hydrogen differential 
scattering cross sections for electron Bessel beams of differ¬ 
ent OAM and transverse momentum, denoted by the Bessel 
beam opening angle a (in mrad). Similar features as for the 
screened Coulomb scattering amplitude are present [27]. 



FIG. 9. (color online) The scattering cross sections for exci¬ 
tation of the 2p2 final state. Note that the peaked maximum 
is smooth when viewed on the much smaller d scale of e.g. 
Fig. lb. The small maximum in the region where 9 < a, al¬ 
though real, is relatively insignificant when compared to the 
cross section’s maximum value. 
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FIG. 10. (color online) The scattering cross sections for OAM transfer events to the 2p± states are shown for positive incoming 
beam OAM, 1. The negative I results are omitted because they are identical to these with 2p+ and 2p- swapped due to the 
symmetry between both situations. Legend is identical to the one in Fig. 8. 
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its sign. This is expected as everything involved is 
symmetric with respect to rotations around the 2 ; 
axis. Non-symmetrical transitions depend strongly 
on the value of 


• The scattering electron wave transfers 0AM to the 
atomic state, and this transfer is reflected in the 
outgoing wave. 


• Taking the plane-wave limit = 0, fcj, —>■ 0), one 
recovers the usual plane wave scattering amplitudes 
(see e.g. Ref. [44] and Fig. lb). 

The results of Sec. Ill C are also shown in Figs. 8-10 (in 
atomic units). The analytical results in Eqs. (33), (34), 
(37), and (36), and the numerically integrated scattering 
amplitudes using Eq. (38) were confirmed to be equal. 
These results and their physical implications shall now 
be discussed. Remember that the differential cross sec¬ 
tion is the probability of an electron being detected at a 
certain scattering angle, which is the primary interpre¬ 
tation that will be given here. Note that the plots show 
only a line profile; the scattering amplitudes shown here 
are all cylindrically symmetrical (up to a possible vortex 
phase factor), resulting in the scattered electrons forming 
doughnut-shaped intensity profiles. 

The elastic differential scattering cross section, shown 
in Fig. 8 , is very similar to that of the screened Coulomb 
potential previously treated [27] . For a non-vortex beam, 
£ = 0 , higher transverse momentum shifts this peak off- 
center, although even for an incoming hollow beam (with 
no low-A:j^ components, generated by e.g. an annular 
aperture), a significant amount of intensity is still scat¬ 
tered on-axis. This can be deduced from the additive 
contributions of various transverse momentum compo¬ 
nents to the scattering amplitude (shown in the topmost 
plot of Fig. 8 ), which additively result in a non-zero on- 
axis differential scattering cross section. For beams con¬ 
taining a phase vortex, i.e. £ 0 , there is no on-axis 

intensity and the angle at which the cross section peaks 
shifts outwards with increasing incoming transverse mo¬ 
mentum. This feature could be exploited as a sensitive 
detector for non 0AM preserving transitions in future 
experiments. Another spherically symmetric transition 
is Is —2s. Its inelastic differential scattering cross sec¬ 
tion (given explicitly by Eq. (34)) has the same shape and 
behavior as the Is —> Is cross section, and is therefore 
not explicitly shown. 

The differential cross sections for ground state excita¬ 
tions to various 2 p substates (quantized along the beam 
direction) are shown in Figs. 9 and 10. The 2 p 2 cross sec¬ 
tion shows a narrow peak that moves away from the beam 
axis, giving way for a zero that appears at the threshold 
transverse momentum defined by the condition Qz = 0 
(which is a prefactor in Eq. (37) and is located at: 


00 = cos 



< a. 


(45) 


When an incoming vortex beam transfers 0AM to an 
atomic electron, the associated differential scattering am¬ 


plitude reflects this, as shown in Fig. 10. In general, 
the lower the outgoing 0AM is in magnitude ]£'], the 
higher the scattering cross section is for low angles. An¬ 
other general feature of these results is that transitions 
where 0AM is given to the atomic electron, ]£'] < ]£], 
have a relatively larger differential cross section for scat¬ 
tering angles 6 < a. Note that not all 0AM needs to 
be transfered for this to be visible, and these differences 
are orders of magnitude larger for higher outgoing 0AM. 
Transitions where 0AM is taken from the atomic elec¬ 
tron, ]£'] > ]£], have larger differential cross sections for 
larger angles, 0 > a. This can be clearly seen by com¬ 
paring the plots in the left column of Fig 10 to those in 
the right one (the top row only shows the latter case, 
]£'] > ]£]). Perhaps the most evident form of this is that 
an £ = 1 beam exciting a Is —>■ 2 p_|_ transition will scat¬ 
ter most electrons on the beam axis and even to 0 = 0 . 
This is possible because for this specific transition, the 
outgoing beam has lost its 0AM, and does not suffer 
from the phase singularity previously forcing its ampli¬ 
tude to zero there. The same happens for an £ = — 1 
beam exciting a Is —>■ 2p_ transition, as the scatter¬ 
ing amplitudes are identical. For higher order incoming 
beams, this change is less dramatic, as the outgoing vor¬ 
tex phase still forces the central scattering to zero. Al¬ 
though it is not displayed clearly on these figures’ scales, 
the (red, long-dashed) plane wave differential cross sec¬ 
tion has the same shape as in Fig. lb, and its central zero 
gets pushed outwards for higher transverse momentum. 
The sharp dip at around a = 0 is exactly the zero of 
fiJ^ 2 p± ill Fig. lb at 0 = 0 , but shifted outwards due 
to the non-zero incoming transverse momentum. It coin¬ 
cides with the peak for the 2p^ differential cross section 
at the same scattering angle. These shifts of intensity can 
be understood intuitively by considering the transverse 
profiles of a vortex beam, where higher 0AM is generally 
paired with a larger spatial extent of the wave function. 
Finally, the large angle scattering of both situations also 
shows a quantitative offset which could also be used to 
differentiate the events. 


V. CONCLUSIONS 

We extended inelastic quantum scattering theory to 
non-trivial incoming electron waves, including orbital an¬ 
gular and transverse momentum. A quick review of the 
textbook theory illuminated some often forgotten, but 
important facts about the hnal state and momentum 
transfer. Two different methods were then applied to 
obtain the vortex scattering amplitudes of inelastic tran¬ 
sitions for a hydrogen-like system. Special attention was 
given to the atomic state’s 0AM and the consequences 
of it being changed by a scattering electron. 

The first method involves the Bessel addition theo¬ 
rem, which leads to unwieldy analytical expressions for 
all scattering amplitudes. Nonetheless, these calculations 
led to simple selection rules when 0AM transfer is in- 
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volved and allows to estimate the regime of scattering 
angles for which they are valid. Additionally, a form of 
0AM reciprocity was shown to exist, tied to the 0AM 
transfer and the central scattering amplitude. 

The second approach, using the Fourier representa¬ 
tion of the Bessel beam, resulted in a purely analytical 
method to obtain hydrogen scattering amplitudes, along 
with a less cumbersome numerical solution using an inter¬ 
mediate result. The scattering amplitudes are influenced 
strongly in the presence of 0AM transfer, even outside of 
the ideal selection rule validity regime (0 = 0). Combined 
with energy filtering, the predicted asymmetry could pro¬ 
vide a means to better separate scattering contributions 
for various final states with distinct 0AM, leading to an 
improved measurement of the final state density. This, 
in turn, would provide atomic-resolution magnetic infor¬ 
mation. Specific Am transitions can be filtered from the 
total scattered intensity using any combination of pre- or 
post-selection of 0AM of the scattering electron, and in¬ 
coming and outgoing transverse momentum (i.e. by lim¬ 
iting collection and convergence angles). The selection on 
transverse momentum is routinely done by choosing de¬ 
tector geometry, selected area apertures, and objective 
apertures. Post-selection on 0AM has yet to be prac¬ 
tically implemented efficiently in electron microscopes. 
Once this is in place, though, the use of 0AM selec¬ 
tors will improve selectivity to certain transitions signifi¬ 
cantly, as a large amount of background signal caused by 
electrons with the “wrong” 0AM will be eradicated. 
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Appendix A: Elastic vortex coulomb scattering 
amplitude 

In the cylindrical expansion of the inelastic scatter¬ 
ing amplitudes, the elastic scattering amplitude appears. 
This was calculated in ref. [27] , and the result is repeated 
here: 


fYi(k';k,i)= - 


2meVo re'' 


n - r2 


rir2 VO + ?’2 


Id 


rl = ql+^^ + (k^-k[)^ 


z z , z , 

O — M + 


Here, Vq is either or 


{k± + k\_) . 

depending on the context. 


and the treatment in this paper deals with unscreened 
Coulomb potentials only, so /j, = 0 everywhere. 
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